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Two-dimensional rotational and vibrational temperature measurements were made at the nozzle 
exit of a free-piston shock tunnel using Planar Laser-Induced Fluorescence. The Mach 7 flow consisted 
predominantly of nitrogen with a trace quantity of nitric oxide. Nitric oxide was employed as the 
probe species and was excited at 225 nm. Non- uniformities in the distribution of nitric oxide in the 
test gas were observed, and concluded to be due to contamination of the test gas by driver gas or cold 
test gas. The nozzle-exit rotational temperature was measured and is in reasonable agreement with 
computational modelling. Non-linearities in the detection system were responsible for systematic 
errors in the measurements. The vibrational temperature was measured to be constant with distance 
from the nozzle exit, indicating it had frozen during the nozzle expansion. 


Introduction 

C RUCIAL to the design of future aerospace vehicles is 
the accurate prediction of physical and chemical effects 
influencing the aerodynamics at hypersonic velocities. Re- 
alistic evaluation of hypersonic phenomena on a particular 
vehicle design can be managed in one of three ways: flight 
tests of actual vehicles; simulation of flight conditions on 
scale models in ground-based test facilities; or numerical 
simulations implementing models of physical and chemical 
processes — known as computational fluid dynamics (CFD). 
The first case is prohibitively expensive and reserved for 
the prototype design. Individual ground-based facilities can 
only simulate flow conditions for a portion of the flight tra- 
jectory. The flows produced by such facilities only partially 
simulate flight flow conditions. The latter method, CFD, is 
a promising alternative with the advent of low-cost, high- 
power computing resources. However, CFD codes require 
validation against experimental data before they can be con- 
fidently used to predict inflight behaviour and aid in vehicle 
design and performance evaluation. The experimental data 
must come from the first two methods mentioned above, 
therefore emphasising the close relationship that exists be- 
tween the alternative approaches. 

The shock tunnel is a ground-based test facility capable 
of producing very high enthalpies and shock- layer tempera- 
tures. It consists of a converging-diverging nozzle attached 
to the end of a shock tube. The performance of the shock 
tunnel can be enhanced by employing the free-piston driver 
technique, 1 thus making it capable of producing chemically 
reacting flow phenomena (real-gas effects). Shock-tunnel 
flows generally only last for a few milliseconds due to the 
impulse nature of the shock tube, and as necessitated by the 
extreme gas temperatures attained. 

Many optical techniques have been used for flow visuali- 
sation in shock tunnels and to complement surface pressure 
and heat-transfer measurements. Schlieren photography, 
shadowgraph and interferometry have been used extensively 
to provide two-dimensional measurements with high tempo- 
ral resolution. However, they can only provide an integrated 
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measurement along the line-of-sight through the flow which 
masks any three-dimensional behaviour in the flow. Pla- 
nar laser-induced fluorescence (PLIF) has been widely used 
for flow visualisation in subsonic and supersonic flows and 
as a combustion diagnostic. 2,3 It provides high spatial and 
temporal resolution, two-dimensional quantitative measure- 
ments and is chemical species specific. PLIF uses a thin 
sheet of laser light to excite a radiative transition in a 
particular chemical species in the flow. The laser sheet inter- 
rogates a thin (<1 mm) cross section through the flowfield, 
allowing PLIF to resolve three-dimensional flow features. 
The short duration of the laser pulse and resulting fluores- 
cence makes PLIF particularly useful for impulse facilities 
where the flow may only last for a millisecond. 

The high-enthalpy flows produced by free-piston shock 
tunnels have a unique set of experimental problems not 
encountered in other flow facilities. Flows generated by 
free-piston shock tunnels have large pressure and temper- 
ature variations which cause large fluorescence signal vari- 
ations due to high collisional quenching rates and changes 
to the spectral overlap integral between the laser and the 
absorption transition. 4,5 Other problems include excessive 
laser-beam absorption which can lead to hole burning of the 
laser spectral profile, 6 spectral interference from other flow 
species and fluorescence trapping. A more severe problem is 
that metallic contaminants which are eroded from the shock 
tube walls produce strong emission in high-temperature re- 
gions in the flowfield. 7 This background emission can easily 
overwhelm the fluorescence signals, so careful spectral and 
temporal filtering of the signal is necessary. This remains 
the greatest obstacle for the application of PLIF to shock 
tunnels. 

In order to validate complex flowfields calculated using 
CFD codes with experimental data from ground-based facili- 
ties, one must have an accurate knowledge of the freestream 
conditions produced by the facility. Measurements of ro- 
tational and vibrational temperatures, static pressure, ve- 
locity, and chemical species number densities are neces- 
sary to completely specify the thermodynamic state of the 
freestream. In the past, CFD codes have been used to pre- 
dict these freestream parameters based on calibrations with 
pitot-pressure data, but little work has been performed on 
testing the validity of these nozzle-flow calculations. In the 
current study, we attempt to address this research shortfall 
by concentrating our efforts on measuring rotational and 
vibrational temperatures for nitric oxide (NO) using PLIF. 

The process of validating a CFD code can be a particu- 
larly arduous one. The methodology and strategies used are 
therefore very important. The process of validation should 
commence with the study of simple flows {i.e., continuum, 
no real-gas effects) and gradually increase complexity of 
the flow {i.e., increase Mach number, add chemical effects). 
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Experiments should be designed to test a single flow phe- 
nomenon while minimizing competing effects that may com- 
plicate data interpretation (cause-and-effect tests). Simple 
models can be used to test specific flow phenomena and in 
some cases this is more desirable. Following this method- 
ology, the flow employed in the experiments presented here 
was a low-enthalpy nitrogen flow. The effects of chemistry 
are removed, simplifying the problem considerably. The vi- 
brational temperature freezes during the nozzle expansion, 
and the temperature at which this freezing occurs has a sig- 
nificant influence on the nozzle-exit rotational temperature. 

For example, consider the current flow conditions cal- 
culated with a quasi-one-dimensional inviscid nozzle-flow 
code (described later in this article) where the vibrational- 
freezing temperature and chemistry can be varied. If we 
calculate the two extreme cases for vibrational freezing, one 
where the vibrational temperature is held in equilibrium 
with the translational temperature throughout the nozzle 
expansion and the other where it is frozen at the nozzle- 
reservoir temperature, we get 560 K and 345 K, respectively, 
for the nozzle-exit rotational temperature. By comparison, 
the effect of switching chemistry between equilibrium and 
non-equilibrium produces only a 20-K difference at these 
operating conditions. Hence, by measuring the rotational 
temperature one can get a measure of how much energy is 
stored in the vibrational modes of the test gas, and therefore 
infer the vibrational-freezing temperature of N 2 . 

Ideally one would like to measure the vibrational temper- 
ature, in this case of N 2 since it is the predominant test 
gas species (for example, using broadband CARS 8 ). Unfor- 
tunately in the current work we can only measure the NO 
vibrational temperature with PLIF. Since NO is a minority 
species in the flow (1%) its vibrational-freezing temperature 
has a negligible effect on the nozzle-exit conditions. It can, 
however, be used to test the NO vibrational-energy transfer 
rates used in the nozzle-flow code. 

In order to use NO PLIF to measure temperature, a small 
amount of NO must be present in the flow. Nitric oxide 
is stable and can be easily seeded into the test gas, 9,10 
however, this is not necessary here because of the high tem- 
perature (4430 K) achieved in the nozzle-reservoir region. 
A small amount of 0* 2 added to the test gas is sufficient to 
produce NO at the nozzle exit without the need for han- 
dling toxic NO. The 0* 2 dissociates in the nozzle-reservoir 
region and recombines with nitrogen atoms to form ~1% 
NO. This value is constant downstream of the nozzle throat 
as the flow becomes chemically frozen. This approach also 
provides an opportunity to test the PLIF technique without 
complicating effects from other flow species ( e.g ., 0 2 ), or 
from excessive beam attenuation caused by high NO con- 
centrations which are typical of air flows. 

A recent PLIF study 11 on an arc-heated facility at similar 
flow conditions to the current experiment produced good 
results, however, due to the running time of several hours, 
the experimenters were afforded the luxury of scanning their 
laser across many excitation lines to build up a spectrum. 
They observed vibrational nonequilibrium which was also 
confirmed by N 2 -CARS experiments. 12 Very little work has 
been performed on impulse facilities due to the difficulties 
it entails. The objective of the current work is to explore 
these difficulties and develop PLIF to the stage where it can 
be used alongside other diagnostics regularly employed on 
shock tunnels. 

PLIF Theoretical Considerations 

Laser-Induced Fluorescence Theory 

The theory of PLIF is well developed 2,3 and only the 
points relevant to the current work will be discussed here. In 
particular, we will present a general theory and then make 
approximations which are relevant for the current applica- 
tion. 

In a typical PLIF experiment, a laser is tuned to an 
optically- allowed electronic resonance of a particular molec- 


ular species. The molecules in a particular electronic- 
rotational- vibrational level mv" J" (typically the ground 
state) are excited to an electronic-rotational-vibrational 
level nv' J ' . Molecules in the excited level nv' J' can then 
emit radiation (fluorescence) which is collected by a detec- 
tor (typically a two-dimensional CCD array). Molecules 
may also undergo collisions with other molecules or atoms 
which cause a transfer of population out of the laser cou- 
pled state nv' J’ to nearby rotational levels. This is known 
as rotational-energy transfer (RET). These collisionally- 
excited nearby rotational levels also fluoresce, resulting 
in a broadband fluorescence signal. Similarly, there may 
also be vibrational-energy transfer (VET) from the ex- 
cited vibrational level. Other collisions may cause some 
of the molecules in the excited state to be de-excited non- 
radiatively (collisional quenching) which produces a corre- 
sponding decrease in the total amount of fluorescence signal. 

The number of fluorescence photons reaching a single 
pixel of the CCD detector is given by 3 


N p = NTfBBj.j.EGQ^e , 

47T 


(i) 


where Nt [cm -3 ] is the number density of absorbers; /b is 
the Boltzmann fraction; Bjijn [s -1 (W/cm 2 /cm -1 ) -1 ] is the 
Einstein absorption coefficient for this particular rotational- 
vibrational transition; E [J] is the energy of the laser pulse; 
G [1/cm -1 ] is the spectral overlap integral; <f> is the fluores- 
cence yield; ft is the solid angle subtended by the detector; 
j] is the detection efficiency; and £ [cm] is the length of the 
interaction volume in the direction of the beam. The inter- 
action volume is defined as that volume in the flow whose 
fluorescence is collected by a single detector pixel. The 
length and breadth of the volume are determined by the 
dimensions of the pixel, and the depth is determined by the 
thickness of the laser sheet. 

Each of the parameters in Eq. (1) shall now be defined 
in more detail. The number density of absorbers can be 
defined by Nt = (px/kT) where p [Pa] is pressure, x is NO 
mole fraction, T [K] is the kinetic temperature and k is the 
Boltzmann constant. The Boltzmann fraction is given by 


/b 


2 J" + 1 
—z ex P 

^total 


te) eXp (fe Tu) 


( 2 ) 


Here, we have indicated the rotational and vibrational con- 
tributions to /b and made the distinction between the ro- 
tational temperature T ro t and vibrational temperature T v ib- 
F j >1 is the rotationally dependent part of the energy of the 
absorbing state and G v » is the vibrational energy. Z to ta i is 
the total partition function. 

The fluorescence yield <f> represents the ratio of de- 
excitation transitions that produce fluorescence photons to 
the total number of de-excitation transitions. It is given by 


<f> = 


A e ff 

A-total T Q 


(3) 


where A e ff [s -1 ] is the effective spontaneous emission rate 
for the collected fluorescence, A to tai [s -1 ] is the sum of spon- 
taneous emission rates for all possible radiative transitions 
from the excited level, and Q [s -1 ] is the collisional quench- 
ing rate. A e ff represents the way that A to tai is modified 
by spectrally selective elements in the detection chain (e.g., 
filters, camera quantum efficiency). Therefore the effective 
emission rate is given by 13 


-4 eff = V T; 


\v" A-v' v n 


(4) 


where T\ v n is the spectral transmission function for the de- 
tection system at the wavelength A and the summation is 
over all vibrational bands of the lower electronic state. 
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The spectral-overlap integral G [1/cm 1 ] represents the 
degree of overlap between the irradiance lineshape L(v) 
[1/cm -1 ] and the absorption lineshape Y (v) [1/cm -1 ]. It 
is given by 

G = J Y(v)L(v)dv (5) 

where Y (u) and L(u) are both normalised such that 
J L(v)dv = 1 and f Y(i/)di/ = 1. The absorption lineshape 
defined above is derived from the definition of the spectral 
absorption coefficient, which is given by 

k v = S 12 Y(is) , (6) 

where S 12 [cm -2 ] is the line-strength factor for the transi- 
tion. With k u expressed in this manner, S 12 contains only 
terms relating to the transition probability and is indepen- 
dent of lineshape effects. The exact form of the lineshape 
function is determined by various spectral-line broadening 
processes. 

Thermometry 

There are many different approaches to LIF thermome- 
try, 2,3,14 but here we discuss only the method used in the 
current experiments. This method is based on the two-line 
planar-thermometry technique of McMillin 9 and Palmer. 10 
Consider the temperature and rotational-level dependence 
of each quantity in Eq. (1). Hence, we have 

N p oc E N t (T) / b (T, J) B(J) G(T, J) <F(T, J) , (7) 

where, for simplicity B = Bj>jn, T = T rot and J = J " . 
Now, if we make fluorescence measurements by exciting two 
different rotational levels (denoted by numbers 1 and 2), 
and take the ratio of the signals, we obtain 

Signal Ratio = ^4 = , (8) 

/V p 2 E 2 -D‘2 J B2 

where C is a constant and we have assumed the quantities 
G(T, J) and <F(T, J) are independent of the excitation level 
J and hence any temperature dependence cancels out in 
the ratio. Also, we have assumed that the number density 
of NO molecules TVr(T) remains constant between the two 
measurements and also cancels. We will justify these as- 
sumptions below. The laser energy E is not dependent on 
J but usually fluctuates from pulse-to-pulse and so is re- 
tained for completeness. Equation (8) shows that the only 
temperature dependence is through the Boltzmann fraction, 
and this is the basis of most LIF thermometry techniques. 
Substituting for / B from Eq. (2) and assuming that each 
transition is within the same vibrational band (i.e., G v >> is 
constant), we obtain 


N pl _ Ei Hi (2 J, +1) -F J2 ) 

JV p2 E 2 B 2 (2J 2 + 1) eXP ( kT t ot 


(9) 


and, solving for T rot , 

Trot = / Fj2 - Fjl)/fc . • (10) 

E 1 B 1 (2J 1 +1) N p 2 J 

In a typical temperature measurement, one measures the 
fluorescence signal N p and laser pulse energy E for each 
rotational level. Then, using calculated values for Fjn and 
Bjtjn , the temperature can be obtained. This assumes that 
the value of C is known or has been determined empirically. 
When more than two lines are used a Boltzmann plot of 
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J> J" (2 J" + 1) 


(ii) 


versus Fjn gives a straight line with slope 
— l/(fcT ro t). The y - intercept is ln(C*) where 


C* = N T G$$dr]lexp(-G v n /fcT vib )/£totai is assumed 
to be a constant independent of J. 

Vibrational temperature can be measured in an analo- 
gous way. Measuring fluorescence images for two or more 
transitions with different v" and plotting 

ln ( Np 

11 + 1) exp(-Fjn /kT rot ) 



versus G v n gives a straight line with slope — l/(fcT v ib). By 
keeping Fjn as constant as possible, its influence is minor, 
and the vibrational temperature is therefore determined in- 
dependent of the rotational temperature. 

Assumptions and Calibration 

The constant C in Eq. (10) can be written such that all 
the quantities related to the detection system are collected 
into a new constant C . Hence, we have 


_ fi 1771^1 
D 2 7?2 ^2 

(13) 

lf NtiGi&i 

(14) 

n T2 g 2 & 2 


For PLIF thermometry in unsteady or turbulent flows, it 
is necessary to use two lasers and two cameras. The two 
transitions are excited in rapid sequence, thus freezing the 
flow motion, and ensuring that the flow conditions remain 
the same for both measurements (i.e., N T i=Nt 2 )- However, 
for steady or repeatable flows it suffices to use one laser and 
camera and assume a degree of flow reproducibility. For 
pulsed facilities it is ideal to make the measurements on a 
single shot as each tunnel run is not perfectly repeatable. In 
the current work the measurements were made over many 
tunnel runs because it was not possible due to monetary 
constraints to use two laser/camera systems. To reduce the 
influence of shot-to-shot flow variations, five runs were aver- 
aged for each transition, and the shock speed was monitored. 
Tunnel runs where the shock speed varied by more than 1% 
from the average were rejected. If a single detection system 
is used then Qir/i = £l 2 rj 2 , and if we assume perfect flow 
reproducibility V T i=Nt 2 , and therefore C' = 1. 

In Eq. (8), we assumed that the ratios $\/$ 2 and G\jG 2 
were independent of temperature. This is a good assump- 
tion if <F and G are independent of J and reabsorption of 
fluorescence is negligible. 15 For the A 2 E + (v' = 0) band of 
NO, the spontaneous emission rates and collisional quench- 
ing rates are independent of rotational quantum number. 16 
Rotational-level- dependent quenching has been observed in 
OH A 2 E + for low-J and temperatures in the range 300- 
1200 K. According to Paul et al. 17 and based on com- 
putational models, it should also manifest in NOA 2 E + at 
J < 5 when temperatures are less than 300 K. J-dependent 
quenching would cause a significant systematic error in 
low-temperature NO thermometry if low-J lines were used. 
Although low-J lines are used in the current work, system- 
atic errors are avoided because the temperatures measured 
are greater than 300 K. Hence we have A e ffi = A e ff 2 and 
Q\ — Q 2 and therefore <3> is also independent of J. 

The overlap integral G depends on the laser lineshape 
L(u) and the absorption lineshape Y(i/). The absorption 
line shape, Y(i /), is affected by collisional broadening and 
shifting, 18,19 Doppler broadening, and Doppler shifts due to 
the flow velocity. Only the Doppler broadening varies with J 
but for transitions in the same vibrational band, u a i ~ v a 2 , 
and so Y (i/) can be considered independent of J. Provided 
the laser-line profile overlaps the absorption- line profile in 
the same way for both transitions, G i — G 2 . The overlap in- 
tegral will vary throughout the flowfield due to variations in 
pressure and temperature, but will remain independent of J. 
For vibrational measurements the Doppler width increases 
by 8% between 225 and 245 nm which causes the overlap 
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X (nm) 226 237 248 259 272 285 


Fig. 1 The NO A 2 E+ X 2 U electronic transition. The 
fluorescence spectral distribution is shown at the bot- 
tom of the figure. The branching ratios were taken from 
Piper. 20 

integral to increase by 2.2%. This produces measured vi- 
brational temperatures which are systematically high by 
+0.4%. From the above considerations and the fact that 
C' = 1, it follows that C — 1. Provided that the assump- 
tions about the NO spectroscopy are valid, then the main 
concern here is that the laser is tuned to the transition in 
the same manner for both lines. 

Ideally one would prefer to test these assumptions by cali- 
brating the whole system at a range of known temperatures. 
A suitable calibration device operating at similar pressures 
and temperatures to the shock tunnel flow is difficult to find. 
A continuously running device might at least enable Boltz- 
mann plots to be obtained where deviations from a straight 
line may indicate systematic errors due to the measurement 
system. In the current work, with no such device available, 
the system was tested by performing measurements in the 
shock-tunnel test section filled with a static NO/N 2 mixture 
at room temperature. The main difficulty this presented 
was due J- dependent attenuation of the laser beam along 
the 150- mm beam path from the entrance window to the 
imaged region. Only when transitions with similar beam at- 
tenuation were chosen did the measured temperatures agree 
with the ambient temperature. 

The above analysis ignores the effect of laser-mode fluc- 
tuations on the overlap integral. The radiation contains 
longitudinal modes and the amount of radiation in each 
mode changes from pulse-to-pulse. If the homogeneous ab- 
sorption linewidth is much greater than the laser linewidth 
these fluctuations would have little impact. All the modes 
would fall under the absorption profile. For a dye-laser 
cavity of 30 cm, the cavity- mode spacing is 0.017 cm -1 , 
which corresponds to 11 modes excited under the 0.18 cm -1 
FWHM laser profile. It should be noted that the fluc- 
tuations may be reduced by normalising the fluorescence 
images by LIF measurements performed simultaneously in 
a separate reference cell. Therefore each laser pulse can be 
corrected for overlap-integral fluctuations. This approach 
was used elsewhere 21 to reduce the standard deviation in 
temperature measurements by a factor of 3 in compari- 
son with energy corrections alone. In the current work the 
pressures are quite low and the freestream homogeneous ab- 
sorption linewidth is ^0.02 cm -1 , only marginally larger 
than the mode spacing. Hence, mode fluctuations cause 
large variations in the measured shock-tunnel fluorescence 
signals. Also, the fluorescence images were normalised only 
by the total laser energy. From 50 tunnel runs the stan- 
dard deviation of the fluorescence signal was found to be 
~25%. Averaging 5 to 8 images for each transition helps 


Table 1 The calculated freestream conditions at 285 mm 
from the nozzle throat. 7 sa t was determined experimen- 
tally. 


Parameter 

Freestream 

Temperature (K) 

396 

Pressure (kPa) 

4.4 

Doppler FWHM (cm- 1 ) 

0.11 

Collision FWHM (cm -1 ) 

0.02 

Collision shift (cm -1 ) 

-0.006 

Overlap integral G (cm) 

3.94 

Collisional quenching rate Q (s -1 ) 

3.8xl0 6 

RET rate (s -1 ) 

4.0xl0 8 

Fluorescence lifetime (ns) 

114 

Saturation irradiance J sa t (kW/cm 2 ) 

160+22 


to reduce the influence of these fluctuations and leads to an 
estimated uncertainty in the rotational temperature of 4% 
(due to laser- mode fluctuations alone). 

An alternative way to calibrate a LIF image is to use a 
known temperature somewhere in the image. In supersonic 
flows, this calibration point can be obtained from the stag- 
nation point on a body or alternatively from the freestream 
temperature, if either is known. Another technique which 
measures accurately at a single point may also be used, for 
example, single- pulse/single-point broadband CARS ther- 
mometry. 8 

Thermometry Strategy for Current Work 

The format for the experiments was as follows. The light 
from a single narrowband laser system was formed into a 
sheet and used to probe the shock tunnel flow. An in- 
tensified CCD camera collected the fluorescence signal at 
right angles to the sheet. Only one image was obtained per 
tunnel run. To make temperature measurements, different 
absorption transitions were probed on separate tunnel runs. 
Generally eight images were produced for each transition, 
and averaged. The temperature thus determined is an av- 
erage of the flow conditions for that series of shock tunnel 
runs. 

The excitation scheme employed involves excitation of NO 
through the A 2 E + X 2 n ( v " = 0) band at 225 nm (see 
Fig. 1). Fluorescence is then collected from the t/=0 — >• 
v"=2,3 and 4 vibrational bands by spectrally filtering the 
fluorescence with reflective filters. Detection at the exci- 
tation wavelength is avoided to prevent contributions to 
the signal from laser scatter. Furthermore, fluorescence 
transitions ending in the less-populated vibrational bands 
(- v " > 1) are preferred to reduce the influence of fluorescence 
trapping , which occurs when fluorescence photons are reab- 
sorbed by molecules between the inter action- volume and the 
detector. Broadband-fluorescence detection was employed 
to collect the fluorescence from all the rotational levels pop- 
ulated by collisions and not just the laser-coupled level. 2 
This has two advantages. First, collecting all the fluores- 
cence gives higher fluorescence intensities. Second, most 
modern lasers have a high degree of polarisation which can 
make the fluorescence from the laser coupled level highly 
anisotropic. 22 The levels populated by RET have reduced 
anisotropy due to depolarising collisions and thus collect- 
ing the broadband fluorescence significantly reduces these 
effects. 

PLIF Transition Selection 

The freestream conditions were calculated using the 
quasi-one-dimensional nozzle flow code STUBE as described 
in the section below on flow modelling. These conditions 
were then used to calculate the various spectroscopic quanti- 
ties that influence the selection of transitions. For example, 
at the estimated freestream temperature and pressure, the 
absorption linewidths are small (0.11 cm -1 ) and therefore a 
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Fig. 2 The T2 free-piston shock tunnel. 


Table 2 Transitions selected for thermometry. 




Transition 

Frequency 

Fj 

G v 


v" 


J" 

(cm- 1 ) 

(cm- 1 ) 

(cm- 1 ) 

a 

0 

Op 

tl2 

2.5 

44069.416 

73.58 

948.66 

b 

1 

Ri 

2.5 

42229.348 

73.41 

2824.76 

c 

2 

s Rn 

8.5 

40562.179 

71.37 

4672.68 

d 

0 

R 2 

23.5 

44382.109 

1045.85 

948.66 

e 

1 

R-2 

23.5 

42516.527 

1035.33 

2824.76 


Table 3 Summary of flow conditions. All values are cal- 
culated except for the shock speed and reservoir pressure 
which were measured. The nozzle-exit conditions are 
given at 285 mm from the nozzle throat. Uncertainties 
are determined by varying the input parameters to the 
codes. 


Nozzle-reservoir conditions (from ESTC): 


u s Shock speed (km/s) 

2.34 

±0.02 

Po Pressure (MPa) 

27.9 

±0.7 

T 0 Temperature (K) 

4430 

±50 

po Density (kg/m 3 ) 

21.2 

±0.5 

ho Enthalpy (MJ/kg) 

5.58 

±0.08 

Nozzle-exit conditions (from STUBE): 


Poo Static pressure (kPa) 

4.36 

±0.2 

Temperature (K) 

417 

±10 

poo Density (kg/m 3 ) 

0.0353 

±0.08 

Uoo Velocity (km/s) 

3.11 

±0.02 

M frozen Frozen Mach number 

7.74 

±0.02 

Vibrational freezing temperatures: 


Tv 2 (K) 

2150 


Ty 2 (K) 

1380 


T v NO (K) 

670 



large number of lines were well isolated. A summary of the 
freestream conditions appears in Table 1. 

Lines were chosen on the basis of: a large value of sep- 
aration in ground state energies to maximise temperature 
sensitivity; isolation from nearby lines; minimal attenuation 
of the laser beam; minimal saturation of the transition; and 
a signal strength above the minimum detectable signal level 
of the detector. The transitions selected are shown in Ta- 
ble 2. All the transitions have very low beam attenuation 
and are well isolated. Transitions with similar Fj values 
were chosen so that vibrational temperature measurements 
could also be made. 

The saturation irradiance J sa t is defined as the laser ir- 
radiance I at which the LIF signal 5 sat has reduced to 
half the value it would have in the absence of satura- 
tion, or its linear value SWear- 3 This can be written as 
S s& t = linear (1/(1 + ///sat ) ) . The influence of saturation is 
minimized by choosing transitions in weak branches and by 


varying the laser energy to achieve the same level of satura- 
tion on each transition (J/J sat ~ 2%). Many lines in the NO 
spectrum consist of overlapping line pairs, where two tran- 
sitions originate from the same lower level but terminate on 
closely spaced upper levels. The separation of these upper 
levels is determined by the spin-rotation splitting constant 
of the A 2 E + state. These line-pairs can be used for ther- 
mometry if saturation is negligible, as in Ref. 10, but in the 
current work they were avoided because negligible satura- 
tion could not be assured. 

Flow Modelling and the Facility 

The nozzle-reservoir conditions were calculated using the 
one-dimensional Equilibrium Shock Tube Code (ESTC). 23 
ESTC uses the initial shock tube fill pressure and tem- 
perature, test gas composition, measured shock speed, and 
measured reservoir pressure to determine the reservoir tem- 
perature. This temperature is then used in the quasi- 
one-dimensional inviscid nozzle-flow code STUBE. 24 One- 
dimensional codes are very useful as engineering tools since 
they provide reasonably accurate solutions in only a few 
seconds. 

STUBE was modified to include a simple vibrational- 
freezing model. The vibrational-relaxation time was cal- 
culated using the theory of Landau and Teller and the 
constants from Ref. 25. Using this relaxation time it was 
possible to determine the point during the nozzle expansion 
at which the vibrational temperature freezes. A sudden 
freezing criterion is imposed and the vibrational temper- 
ature held constant for the remainder of the expansion. 
The model accounts only for vibrational-translation (V- 
T) energy transfers and vibrational-vibrational (V-V) en- 
ergy transfers are ignored. Due to limited availability of 
relaxation-rate data, the Landau- Teller rate constant for 
NO-N 2 collisions was assumed to be the same as that for 
NO-NO collisions. However, the relaxation rate for NO-NO 
collisions is greater than that for the NO-N 2 collisions. 26 
The model therefore overestimates the relaxation rate of 
NO and thus under predicts the NO vibrational freezing 
temperature. This has little effect on the other calculated 
nozzle-exit parameters since NO is a minority species in the 
current experiment. It does, however, provide a point of 
comparison with the experimentally measured NO vibra- 
tional temperature. 

Figure 2 shows a schematic of the T2 shock tunnel. 1 A 
1.20- kg piston is free to move inside the 3.1-m compres- 
sion tube which has an internal diameter of 76 mm. The 
shock tube is 1.98 m in length and has an internal diam- 
eter of 21 mm. Attached to the end of the shock tube is 
a converging-diverging nozzle which exits into a test sec- 
tion which has optical access through 4 window ports. For 
the current experiments, the primary diaphragm was made 
from 0.74-mm-thick mild steel which has a burst pressure 
of 46.9 MPa. The secondary diaphragm was made from 
0.025-mm-thick mylar. The present study employed a 7.5° 
half-angle axisymmetric conical nozzle. The throat of the 
nozzle is 6.35 mm in diameter, the axial length of the ex- 
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Fig. 3 Experimental apparatus for LIF imaging experiments. 


panding section is 255 mm and the nozzle-exit diameter is 
73.6 mm. A pitot-pressure survey was conducted at the noz- 
zle exit to estimate the size of the inviscid core. A reduced 
nozzle angle of 7.0° was used with STUBE to account for 
the effect of the boundary layers on the nozzle walls. 

The tunnel was operated with a helium/argon driver 
(81.3 kPa He, 37.3 kPa Ar) to obtain tuned- piston oper- 
ation. 27,28 The shock tube was filled with 100 kPa of a 
mixture of 1.1% O 2 in N 2 at 293 K. The piston is driven 
by air from the high-pressure reservoir which is filled to 

4.34 MPa just prior to the shot. When the piston is released 
it is propelled down the compression tube, compressing the 
driver gas and causing the diaphragm to rupture. A shock 
wave propagates through the test gas in the shock tube and 
its speed is measured by three timing transducers along its 
length. The transducers are placed at 1281, 316 and 11 mm 
from the nozzle inlet. The speed computed based on the 
shock transit time between the first two transducers was 
2.49 km/s and between the last two it was attenuated to 

2.34 km/s due to viscous effects. The shock speed is ex- 
pected to decrease rapidly near the start of the shock tube 
and be relatively constant by the time it gets to the nozzle 
inlet. Hence, the value of 2.34 km/s was used in calculating 
the nozzle-reservoir conditions. The shock wave reflects at 
the shock tube end producing a high pressure and temper- 
ature reservoir of test gas. The measured nozzle-reservoir 
pressure was 27.9 MPa. Prom these values the reservoir tem- 
perature was calculated using ESTC to be 4430 K, and the 
total enthalpy was 5.8 MJ/kg. The conditions calculated 
using STUBE at 285 mm from the nozzle throat are shown 
in Table 3. 

Experiment 

The experimental arrangement is shown in Fig. 3. The 
laser source is a frequency- doubled excimer-pumped dye 
laser producing up to 6 mj of radiation at 225 nm, with 
a linewidth of 0.18 cm -1 . The pump laser was a XeCl ex- 
cimer laser (Lambda Physik, EMG 150 ETS) operated as 
an unstable resonator, and produced approximately 250 mj 
at 308 nm. The dye laser (Lambda Physik, Scanmate II) 


was operated with different Coumarin laser dyes (Lamb- 
dachrome LC4500, 4700, 4800) to produce 40 mj of tunable 
narrowband radiation between 450 and 490 nm. This was 
then frequency doubled in a BBO I crystal to produce about 
6 mj of UV light between 225 and 245 nm. Approximately 
10% of the laser radiation was split off with a beamsplitter 
(uncoated fused silica) and directed through a small tur- 
bulent flame to perform LIF excitation scans. The NO is 
produced by entrainment of the air in the flame. The LIF 
was collected at right angles to the beam and then imaged 
onto a 0.5-m spectrometer. The spectrometer was operated 
as a 9-nm bandpass filter which allowed detection of LIF 
from the NO (0,1) or (0,2) vibrational bands while blocking 
the very strong OH fluorescence at 308 nm. The flame was 
chosen (rather than a reference cell containing NO) so that 
excitation scans of the v"=l and 2 vibrational bands of NO 
could also be performed. 

The remaining laser radiation passes through a periscope 
to bring the beam to the same level as the shock-tunnel test 
section. An 80-mm wide sheet was produced using a 20- 
mm-focal-length cylindrical lens and a large spherical lens 
(1000-mm focal length, 100- mm diameter). The variation in 
energy across the laser sheet was measured simultaneously 
by use of a dye cell (100-mm length) and a standard CCD 
camera (Pulnix TM-760). 29 This energy distribution mea- 
surement was used to correct the PLIF image obtained on 
the same laser pulse. The laser sheet was focussed in front 
of the test section so that the sheet thickness at the test sec- 
tion was 0.8 mm. This reduced the possibility of saturation. 
The total beam energy was monitored with a back reflection 
from a prism using a UV-sensitive photodiode. The cylin- 
drical lens overexpands the beam considerably so that the 
total sheet energy was reduced to only ~60 /iJ. This corre- 
sponds to an irradiance of 3.7 kW/em 2 , where we have used 
a sheet area of 0.64 cm 2 and assumed a laser pulse duration 
of 25 ns. 

The saturation irradiance was measured using a static 
NO/N 2 mixture in the shock tunnel test section, as de- 
scribed in the section on calibration. Making allowance for 
the change in conditions between this static measurement 
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Fig. 4 Schematic of nozzle exit. The box indicates 
the region imaged by the camera. The imaged region 
is 50x52 mm. The image shown is a rotational tem- 
perature map which has been reflected about the flow 
centreline. 

and the freestream conditions, the estimated saturation ir- 
radiance for the freestream is 160 kW/cm 2 . Therefore the 
laser irradiance used for the experiment was only 2% or the 
saturation irradiance (I /la at ~ 0.02). 

The fluorescence from the test section was collected at 
right angles to the sheet and reflected from two mirrors 
and into the intensified CCD camera. The mirrors are 
reflective filters, designed to maximise the fluorescence sig- 
nal and minimize flow luminosity. The detector used for 
these experiments was an in-house assembled intensified 
CCD camera. 30 The CCD camera (EEV Photon, P/N 
P45580/V5.3/PHO) is an 8- bit asynchronous device which 
is triggered separately from the intensifier. The fluorescence 
lifetime in the freestream is ~-T14 ns and so an intensifier 
gate of 650 ns was used. This allows greater than 99% of 
the signal to be collected and about 40 ns between the start 
of the gate and the laser pulse to account for time jitter 
in the triggering system. Several ‘laser off’ shots were also 
performed to verify that there was no flow luminosity. 

Procedure 

To ensure an accurate measurement of the laser energy 
used in the experiment, it was necessary to ensure that 
the test-section windows were cleaned prior to each tun- 
nel run. The tunnel was then reloaded and pumped to less 
than 25 Pa. A LIF excitation scan was then performed in 
the flame which ensured that the laser was tuned to the 
correct NO transition. The laser was run at 1 Hz while the 
tunnel was being filled. Immediately before firing (< 2 s), 
the tunnel operator stopped the laser via a remote switch 
next to the firing valve. After the firing valve was opened, 
the nozzle-reservoir pressure transducer detected the shock 
arrival at the end of the shock tube and the laser was fired 
350 jjs later. This delay was chosen to coincide with flow 
establishment in the nozzle and prior to the expected ar- 
rival time of driver gas. The intensified camera and the 
sheet monitoring camera were also triggered from the nozzle- 
reservoir pressure transducer after a suitable delay. 

Results and Discussion 

The imaged region is defined by the box in Fig. 4 which 
shows a rotational temperature map at the nozzle exit. The 
flow is from left to right. For these experiments, only half of 
the freestream was imaged, and this 50x52- mm region in- 
corporates the flow centreline. The camera actually views a 
55-mm wide region, but the laser sheet is reduced to 50 mm 



Fig. 5 A collection of LIF images showing flow non- 
uniformities. The position of each image relative to the 
nozzle is shown on Fig 4. 



Fig. 6 Averaged LIF images of the freestream. The 
flow is from left to right. Images (a), (b) and (c) are for 
the low-J lines and show the same characteristic shape. 
Images (d) and (e) are for the high-J lines. Image (f) 
is the signal produced when the laser is detuned from a 
transition at 245 nm. Image (c) includes the correction 
for the non-resonant background (f). 

so that the edges of the sheet observed in the LIF image can 
be correlated with the sheet-energy correction image. The 
imaged region starts at 16 d= 2 mm from the nozzle exit and 
includes a portion of the centreline as well as the edge of the 
core flow. Here, allowance has been made for the recoil of 
the shock tunnel when it is fired, which causes the nozzle to 
move 8 ± 1 mm further away from the imaged region. Note 
that the bottom part of the image outside the box is simply 
the reflection of the top portion of the image. 

For a uniform freestream flow, a uniform fluorescence 
signal is expected. However, during the course of the ex- 
periments, several flow non- uniformities were observed (see 
Fig. 5), particularly along the centreline of the flow. Regions 
of no LIF signal were observed in approximately one third 
of the 68 tunnel runs performed. They have the appearance 
of small turbulent clouds of gas. 

Several explanations exist for these nonuniformities. 
First, these dark regions may be explained as pockets of 
driver gas that have arrived prematurely, perhaps due to 
instabilities at the contact surface or jetting of driver gas 
into the nozzle-reservoir region via shock wave/boundary 
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Distance from centreline (mm) 

Fig. 7 Cross sections through the rotational tempera- 
ture map perpendicular to the flow direction. 

layer interaction. 31,32 The driver gas has no NO in it and 
hence produces no LIF signal. Second, because the NO 
used in the current experiment is generated in the shock- 
reflection process, there exists a boundary layer of test gas 
on the shock-tube walls without NO. This boundary layer 
may mix with the test gas in a similar method to the shock 
wave/boundary layer interaction described above. Third, 
the lack of signal may be due to particulate matter obstruct- 
ing the fluorescence between the plane of the laser sheet and 
the exit window. It is known that the nozzle-boundary layer 
entrains particulate matter from the nozzle walls since this 
has been observed by Mie scatter when detecting at the laser 
wavelength (see Fig. 6(f)). However, if this was the cause 
then one would also observe sudden drop outs in the laser 
beam which would cause streaks in the fluorescence signal. 
These are never observed so the it is more likely one of the 
first two scenarios is correct. Attempts to seed the driver gas 
with NO produced negative results which are inconclusive. 

It is interesting to note that very little mixing has oc- 
curred, as indicated by the high contrast of some of the 
‘dark-gas’ pockets. At later test times (>350 /is), LIF im- 
ages display a gradual uniform decrease in intensity which 
may be due to either a drop in the nozzle-reservoir pressure 
and hence nozzle-exit density, or an increased proportion of 
driver gas uniformly mixed with the test gas, or both. 

Individual images showing flow non-uniformities were re- 
jected from the temperature analysis. Each image obtained 
from a single tunnel run was corrected for camera back- 
ground offset, absolute laser energy, and the variation of the 
laser energy across the image. The corrected images for each 
transition were then averaged, typically 8 images per tran- 
sition. These averaged images are shown Fig. 6. Because 
of the limited dynamic range of the camera and the signif- 
icant variation in fluorescence signal intensities, especially 
from different vibrational bands, the camera intensifier gain 
was varied for different transitions. The intensifier gain con- 
trol was therefore calibrated at the same time that linearity 
tests were performed on the intensifier. This explains the 
intensity-scaling factors at the bottom right of each image 
in Fig. 6. For the six images (a-f) in Fig. 6 the gains were 
6.5, 7.2, 7.5, 6.8, 7.2 and 7.5, respectively. Images taken at 
higher gains, such as image (c), show significant shot noise 
despite averaging. The reflective filters efficiently block the 
laser wavelength when excitation is performed in the (0,0) 
and (0,1) bands. However, elastically- scattered laser light 
is not filtered out by the detection system for excitation of 
the (0,2) band near 245 nm. To account for this, several 
runs were performed with the laser at 245 nm but detuned 
from any NO transitions. The average of 6 such images is 
shown in Fig. 6(f). The signal occurs predominantly at the 
edge of the core flow and is most likely due to Mie scatter- 
ing from particulate matter that has been entrained in the 
nozzle- wall boundary layer. There is no signal in the central 



Fig. 8 Rotational temperatures along the centreline of 
the flow compared with STUBE calculations. 

core of the flow indicating it is relatively clean of particu- 
late matter. Image (c) is a corrected LIF image obtained by 
subtracting image (f) to account for the elastically-scattered 
laser light. 

A rotational temperature map was calculated by taking 
the ratio of images (a) and (d) in Fig. 6 and using Eq. 10. 
The temperature map produced is shown in Fig. 4. However, 
a more quantitative assessment can be made by examining 
the vertical cross sections given in Fig. 7. Each cross sec- 
tion was made by averaging a 50-pixel-wide vertical slice 
through the image. The pixel density was 9.2 pixels/mm, 
and so this corresponds to approximately a 5- mm wide slice. 
The data was further smoothed by running a 20- pixel- wide 
(2 mm) integration window through the data. The cross 
sections are at 10-mm intervals and clearly show the de- 
crease in temperature with distance from the nozzle exit. 
Note that the data shown in Fig. 7 greater than +10 mm 
from the centreline is a reflection of the data from below 
-10 mm. The first slice, corresponding to 274-mm from the 
nozzle throat, shows a marked temperature rise at the edge 
of the core flow. This temperature rise is caused by vis- 
cous dissipation 33 in the nozzle- wall boundary layer, and on 
closer examination of the image in Fig. 4, can be seen to 
agree very well with the position of the nozzle wall. There 
is a small temperature rise on the centreline which may be 
caused to some upstream disturbance. Slices taken further 
downstream show a more uniform temperature profile. A 
second temperature map was made from transitions in the 
(0,1) band (images (b) and (e) of Fig. 6) and showed similar 
characteristics, however the signal-to-noise ratio was signif- 
icantly poorer. 

Figure 8 shows streamwise temperature measurements 
compared with STUBE calculations. Each experimental 
data point was obtained by averaging over a 2 x 10 mm 
region at several points along the flow centreline. The two 
limiting cases for vibrational relaxation are shown to indi- 
cate the confidence interval for the calculations. The upper 
curve represents the case when the vibrational temperature 
is equilibrium with the kinetic temperature while the lower 
curve indicates the extremely unlikely case where the vibra- 
tional temperature is frozen at the nozzle-reservoir value. 
The solid curve represents the case obtained by letting the 
N 2 vibrational temperature freeze at 2150 K during the noz- 
zle expansion, as explained above. Good agreement with 
STUBE is obtained with the temperature from the v"=0 
level but the temperature from the v ff =l level is high by 
30 K. 

Temperature Uncertainties 

The major uncertainty in the freestream measurement 
is due to fluctuations in the laser spectral profile. More 
precisely, it is due to the fluctuation in the distribution of 
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Fig. 9 Boltzmann plot for vibrational thermometry. 


energy amongst the cavity modes of the laser. The uncer- 
tainty in a single image due to laser-mode fluctuations was 
determined by considering the data from all tunnel runs to 
be 25%. Averaging 8 images gives an uncertainty of 9.4% 
for each averaged image in Fig. 6. For two-line thermome- 
try, we use the ratio of two such images and therefore the 
uncertainty in the ratio is \[2 greater. The uncertainty in 
the temperature ST is related to the signal ratio R according 
to 


ST_ _ kT SR 
~T ~Ke~R ’ 


(15) 


where A E is the difference between the energy levels (A E = 
Fj i — Fj 2 ). For the freestream rotational-temperature mea- 
surements A E = 972 cm -1 , T 450 K and A E /kT « 3.25, 
and therefore the uncertainty in the rotational temperature 
due to laser- mode fluctuations is only 4%. Other sources of 
random error include the laser-energy measurements, cam- 
era gain and gate repeatability, time jitter in the laser 
system, laser-tuning inaccuracies and flow repeatability. All 
of these uncertainties are negligible in comparison with the 
laser-mode fluctuations. 

Systematic uncertainties have been reduced or avoided 
by careful experiment design and transition selection. We 
estimate that beam attenuation, signal interference and sat- 
uration produced a combined error of less than 1 K. It 
is possible that transition overlaps or other interferences 
may occur that are not predicted by the calculations, but 
these can only be discovered by temperature verification 
experiments. Another systematic error source is due to non- 
linearities in the intensified camera system. This a possible 
explanation for the different rotational temperatures derived 
from the v"=0 and v"=l bands. The linearity of the detec- 
tor was tested by performing LIF in a static N 2 /NO mixture 
and plotting the integrated LIF intensity against laser en- 
ergy. However, the test was ineffective at gains above 7.0 
because shot noise caused a large number of the pixels to 
saturate (exceed the 8- bit dynamic range of 255 counts). 
The test could only indicate the camera was linear up to 
100 counts. Because the signal from the low- J line used for 
the v”=l rotational temperature measurement is stronger 
than the high- J line, it is more likely to suffer from satura- 
tion effects. A reduced low- J signal would produce a higher 
rotational temperature as observed. The signals obtained 
from the v" =0 band were measured at lower gains where 
the camera was verified to be linear up to 220 counts. To 
get agreement between the two measured temperatures the 
signal from image (b) would need to be 20% greater. 

Vibrational Temperatures 

Freestream vibrational temperatures were also deter- 
mined from the LIF images shown in Fig. 6. Using the 
signals from the U'=0,1 and 2 low-J transitions (a), (b) and 
(c) and plotting the expression in Eq. 12 versus G v n gives 


Fig. 10 Vibrational temperatures along the centreline 
of the flow. 

a Boltzmann plot, as shown in Fig. 9. From the slope of 
this plot we can get the vibrational temperature. A similar 
approach is taken to that of the rotational measurements 
to reduce shot-noise uncertainties by averaging the tem- 
perature over a 2 x 10 mm region. Figure 10 shows the 
vibrational temperatures at various distances from the noz- 
zle exit along the flow centreline. The temperature error 
bars are due to laser-mode fluctuations. It is observed from 
Fig. 10 that the measured vibrational temperature is signif- 
icantly higher than that predicted from the computational 
model. The temperature is constant with distance from the 
nozzle exit indicating that it is frozen. The measured tem- 
perature is 785 + 30 K compared with the value from the 
vibrational relaxation model of 670 K. The underprediction 
of the vibrational freezing temperature was expected due to 
the lack of vibrational relaxation data. 

The estimated systematic uncertainties for the vibrational 
temperature measurements due to beam attenuation and 
signal interference from nearby lines are ST /T « +0.7% 
and « +0.3%, respectively. Uncertainties due to satura- 
tion are negligible because the laser energy was varied to 
achieve the same amount of saturation for each transition. 
The Doppler width decreases by 8.0% between 225 and 
245 nm, which causes the overlap integral to increase by 
2.2%. This produces measured temperatures that are sys- 
tematically high by +0.4%. Combining these uncertainties, 
the estimated systematic error for the vibrational measure- 
ments is +0.85%, or +7 K. 

However, it is also believed that camera nonlinearities 
have contributed to lowering the measured vibrational tem- 
perature in a similar way to that described for the rotational 
temperature measurement. From the Boltzmann plot in 
Fig. 9, it can be seen that the signal intensity from the 
v" —2 level is ~35% less than that predicted by a straight 
line fit from the v" =0 and 1 signals. Since this signal was 
acquired at the highest intensifier gain it is most affected 
by camera saturation. A lower v" =2 signal would lead to 
a lower observed temperature. In fact if we use only the 
first two vibrational levels (U'=0,1) we get a temperature of 
^850 K. However, just as with the rotational measurements, 
we suspect that the v ,f =l transitions may also be affected 
by camera saturation, so the actual vibrational temperature 
may be even higher. Increasing the v"=l signal by 20% as 
described in the rotational measurements above leads to a 
temperature of ^950 K. 

Agreement with CFD 

To calculate the nozzle-exit conditions accurate knowl- 
edge of the nozzle-reservoir temperature is required. Here 
we have used the average shock speed measured between 
two transducers 11 and 316 mm from the nozzle inlet. We 
have assumed that this is a good estimate for the shock 
speed at the nozzle inlet. If the shock speed is decreas- 
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ing at the linear rate along the shock tube than the shock 
speed at the nozzle inlet would be only 2.25 km/s, leading 
to an underestimate of of the nozzle-reservoir temperature 
by 200 K, and nozzle-exit rotational temperature by 30 K. 
Furthermore, the calculations for the nozzle-reservoir con- 
ditions (using the code ESTC 23 ) make several assumptions 
and exclude losses due to viscous, radiation and conduction 
effects, driver gas contamination and shock-tunnel impuri- 
ties. 

If the N 2 vibrational-freezing temperature could be mea- 
sured, for example with broadband CARS, then it would 
be possible to make a reasonable estimate for the nozzle- 
reservoir temperature. This would be one method of vali- 
dating the assumptions described above used by ESTC. In 
related work we measured the temperature at the stagna- 
tion point on a hemisphere using PLIF. This work will be 
presented elsewhere, but we mention it here since a mea- 
surement of the stagnation point temperature would give a 
lower bound for the reservoir temperature. The vibrational 
temperature of NO at the stagnation point on a hemisphere 
was measured to be 4500±270 K. Unfortunately the large 
uncertainty in this measurement prevents a definitive es- 
timate for the reservoir temperature, but it demonstrates 
an alternative approach which may aid in validating shock- 
tunnel CFD codes. 

Conclusions 

Planar laser- induced fluorescence of NO has been used to 
characterise the nozzle flow from a free-piston shock tunnel. 
Flow non-uniformities have been observed which indicate a 
contamination of the test gas by driver gas or test gas from 
the boundary layers on the shock tube walls. This is the first 
time to the author’s knowledge that such non-uniformities 
have been observed. This discovery may have significant 
consequences for free-piston shock tunnel research, in par- 
ticular in experiments involving combustion research, where 
a pocket of driver gas passing through a combustor could 
extinguish or considerably affect ignition. 

A rotational temperature image was made at the nozzle 
exit and shows reasonable agreement with CFD calcula- 
tions. Temperatures of 440 ± 25 K and 460 d= 25 K were 
measured at 285 mm from the nozzle throat for transitions 
in the v" =0 and 1 vibrational levels. Non-linearities in the 
homemade camera system are believed responsible for incon- 
sistencies between these two measured temperatures. The 
first measurement was performed at lower intensifier gains 
where the system’s linearity had been validated. The vibra- 
tional temperature of NO at the nozzle exit was measured to 
be 785 ±30 K. It is believed that this temperature may actu- 
ally be as high as 950 K and that camera non-linearities have 
caused a systematic error in the measurement. Nonetheless, 
the temperature was found to be constant with distance 
from the nozzle exit indicating that it is frozen. 

We conclude that PLIF is a useful technique for observing 
flow non- uniformities and mapping flowfield distributions. 
More work is required to validate its accuracy, however 
the problem with camera non-linearities is easily solved 
with modern commercial camera systems. Validation of 
nozzle- flow codes requires improved knowledge of the nozzle- 
reservoir conditions. 
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